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A formal relation between the GT part of the nuclear matrix elements Mgy of Ov/3/3 decay and 
the closure matrix elements M^f of 2u/3j3 decay is established. This relation is based on the integral 
representation of these quantities in terms of their dependence on the distance r between the two 
nucleons undergoing transformation. We also discuss the difficulties in determining the correct 
values of the closure 2u/3/3 decay matrix elements. 
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I. INTRODUCTION 

The 2^/3/3 decay mode of most candidate nuclei has been observed and its half-life determined. Thus, 
the corresponding nuclear matrix elements M 2l> are known. These matrix elements, of dimension energy -1 , 
vary abruptly between nuclei with different Z and A; they exhibit pronounced shell effects. In contrast, 
the fundamentally more important fV/3/3 decay mode have not been reliably observed so far. Therefore, the 
corresponding dimensionless nuclear matrix elements M " must be evaluated theoretically. These calculated 
J> ' quantities, whether based on the QRPA nuclear shell model 043' or the Interacting Boson Model [8|, 

do not show such a variability; instead they vary relatively smoothly between nuclei with different Z , and A. 
Evaluation of the M 0v using the Generator Coordinate Method Q or the Projected Hartree-Fock-Bogolyubov 
Method vary smoothly with Z and A as well. 

However, the calculated values of the matrix elements M Ql/ using different approximations do not agree with 
each other perfectly, differences of a factor of about two exist. Moreover, given the fundamental importance 
of these quantities for the planning and interpreting the Oi//3/3 decay experiments, it would be good to have 
independent observables that could be linked to their magnitude. In that context we wish to address in this 
work several questions: 
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Can one understand intuitively the different behavior of A1 0l/ and M 2u when of Z and A are varied ? 



H • Is there a formal relation between M 0v and M 2v ? 

• If such relation exists can it be used to test the calculated values of the 0^/3/3 matrix elements ? 

In discussing these problems we basically follow the work published by us and our collaborators earlier in 
Ref. [n|. 

II. FORMALISM 

Assuming that the 0uj3/3 decay is caused by exchange of the light Majorana neutrinos, the nuclear matrix 
element consists of Gamow- Teller, Fermi and Tensor parts, 

M°» = M°J T --?- + M°» ee M^ T (1 + xf + XT) , (1) 
9a 
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where \f and \t are matrix element ratios that are smaller than unity and, presumably, less dependent on 
the details of the applied nuclear model. 

In the following we concentrate on the GT part, Mq t , which can be somewhat symbolically written as 

Mgt = (/|E Jfc *j ■ 5 k T+T+H(r lk ,E)\i) , (2) 

where H(ri k ,E) is the neutrino potential and ri k is the relative distance between the two neutrons that are 
transformed in the decay into the two protons. 

The dependence of M 0v on the distance ri k is described by the function C°"(r) (first introduced in see 
also 0,0) 

P OO 

C°J T (r) = (f\E lk a t ■ d k T+T+5{r - r lk )H{r lk ,E)\i) , M° G V T = / C°J T (r)dr , (3) 

Jo 

where S(x) is the Dirac delta function. 

In analogy to Eq. ([3]) we can define for the case of the 2vfij3 decay a new function 

/>oo 

C 2 J(r) = (m k ai ■ a k S(r - ri k )r+r+\i) , M% = / C 2 ?(r)dr . (4) 

Jo 

This function, therefore, is related to the dimensionless closure matrix element for the 2^/3/3 decay, not the 
true, dimension energy" 1 matrix element M 2v that contains the corresponding energy denominators. While 
the matrix elements M 2v and M 2 " get contributions only from the 1 + intermediate states, the function C 2 Y(r) 
gets contributions from all intermediate multipoles. This is the consequence of the 5 function in the definition 
of C^(r). Naturally, when integrated over r only the contributions from the 1 + are nonvanishing. 
It is now clear that, by construction, 

C°J T (r) = H(r,E)xC 2 J(r) , (5) 

which is valid for any shape of the neutrino potential H(r, E). Thus, if Cj'(r) is known, Cq T (r) and therefore 
also Mq^ can be easily determined since the neutrino potential H(r,E) is known and only weakly dependent 
on the average excitation energy E. The equation ([5]) represents the basic relation between the Of and 2v 
/3/3-decay modes. 

Examples of the function C^(r) are shown in Fig. Q]. Note that while the function C^(r) has a substantial 
negative tail past r ~ 2 — 3 fm, these distances contribute very little to Cgy(r) and, therefore also to M 0l/ . 
This is a consequence of the shape of the neutrino potential H (r, E) that decreases fast with increasing values 
of the distance r. 



III. RESULTS AND DISCUSSION 



Remembering that in a nucleus the average distance between nucleons is ~1.2 fm we can somewhat schemat- 
ically separate the range of the variable r in the functions C(V T (r) and C 2 "(r) into the region r < 2-3 fm 
governed by the nucleon-nucleon correlations, while the region r > 2-3 fm is governed by nuclear many-body 
physics. From the form of Cgy(r) it follows that the matrix elements M(V T are almost independent of the 
"nuclear" region of r and hence one does not expect rapid variations of their value when A or Z of the nucleus 
is changed. On the other hand, the 2v closure matrix elements M 2 Y depend sensitively on that region of r 
since there is a substantial cancellation between the positive part at r < 2 — 3 fm and the negative tail at 
r > 2 — 3 fm. Hence one expects sizable shell effects, i.e. a significant variation of M 2v and M 2 f with A and 
Z, in agreement with observations. We have thus answered the first two questions posed in the Introduction. 

But answering the third question is not so simple. While the nuclear matrix elements M 2l> are simply related 
to the 2v half-life T 2 J 2 , and are thus known for the nuclei in which T 2 J 2 has been measured, the closure matrix 

elements M 2 " need be determined separately. One can rely on a nuclear model (e.g. QRPA or nuclear shell 
model) , adjust parameters in such a way that the experimental value of M 2v is correctly reproduced, and use 
the model to evaluate M 2 ". Alternatively, one could use the measured and f3 + strength functions and 
assume coherence (i.e. same signs) among states with noticeable strengths in both channels. In this way an 
upper limit of M 2 " can be obtained. Obviously, neither of these procedures is exact. In the following we leave 
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FIG. 2: Running sums of M^i for different selected nuclei as a function of the excitation energy of the 1 + states in the 
intermediate odd-odd nucleus. 



temporarily aside the important question to which extent the function C^"(r) is strongly constrained by its 
integral value M^". 

We have shown in Ref. [ll[ that, within QRPA, the quantities are negative in many cases, contrary 
to simple expectations (because one obviously expects that the average excitation energy is positive). On the 
other hand, the second method, mentioned above, leads to the positive M^f by definition. Thus, the two 
methods disagree with each other. 
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Since there is a substantial experimental activity devoted to the determination of the /3 strengths [13], it is 
worthwhile to examine in more detail the somewhat unexpected finding that in many cases M 2v and M 2 " have 
opposite signs. Obviously, this has to do with the different weights of the individual 1 + intermediate states in 
M 2v and M 2 Y . Because of this the higher energy excited 1 + states contribute substantially more to M 2 " than 
to M 2v . In Fig. [2] we plot the running sums of M 2 " as the function of excitation energy in the intermediate 
nucleus. One can see that there are negative contributions to the M 2 " values that arise from excitation energies 
E ex > 10 MeV where it is difficult to explore experimentally the corresponding (3 + strength. Note, moreover, 
that this is the region of the giant GT resonance, with substantial fir strength. Hence even a relatively small 
fi + strength can have a sizable effect. 

The negative contribution to M 2v and M 2 f from relatively high lying 1 + states seems to be a generic effect 
present in essentially all nuclei. (Note that analogous effect is present in the shell model evaluation of M 2v for 
48 Ca, see Ref. (l4|. Fig. 1.). Thus, it appears that in many cases when the matrix elements M 2v (with energy 
denominators, see [ll() are evaluated as a function of the intermediate state excitation energy, the final correct 
value is reached twice; once at relatively low excitation energy and the second time asymptotically. 

These negative contributions to M 2v and M 2 f, coming from the vicinity of the giant GT state, clearly exist 
only in calculations that are able to describe the giant GT resonance. In the restricted shell model calculations 
(without the full set of the spin-orbit partners) such high lying 1 + will be absent. Thus, we need to make 
sure that a noticeable fi + strength connecting the final nucleus with the 1 + states in the vicinity of the GT 
resonance really exists and is not an artifact of our QRPA evaluation. Unless and until this dilemma is resolved 
it is premature to proceed further in our original program of finding connection between the nuclear matrix 
elements of the Ov and 2v /3/3-decay modes. 

Nevertheless, we have established a formal relation between the nuclear matrix elements M^ T of the neu- 
trinoless fifi decay and the MB", the closure matrix element for the 2vfifi decay. We also pose a challenge to 
both experimentalists studying the charge exchange reactions of the (n,p) type, and to theorists using methods 
alternative to QRPA to establish whether a noticeable fi + strength at energies near the giant GT resonance 
exists or not. 
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